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Abstract. 

The amplitude of jet distortions and accompanying pressure and velocity 
0\ , fluctuations resulting from Kelvin-Helmholtz instability of three dimensional 

^^ , relativistic jets are explored. The effect of instability on jets as they accelerate from 

VQ ■ sub- to super- Alfvenic speeds is explored and a relativistic stabilization mechanism 

r~^ ' for trans-Alfvenic jets is proposed. The level to which asymmetric instabilities 

^^ ' on supermagnetosonic relativistic jets will grow is predicted theoretically and a 

Doppler boosted "apparent" emissivity is computed. Effects due to helically twisted 
filamentary structure produced by asymmetric modes of instability should be readily 



r^ . observable on relativistic jets. 

o 

C/3 i 1. Introduction 

There is compelling evidence for relativistic highly coUimated jets in both galactic 
(^ ' and extragalactic objects 0. In the galactic "superluminals" GRS 1915-1-105 |^ 

i^j] , and GRO J1655-40 Q l5| the observed proper motions indicate jet flow speeds on the 

C^ ' order of 92% of light speed. Among the extragalactic superluminals, observed motions 

indicate jet flow speeds as large as 99.9% of light speed, e.g., 3C 345 Q. It is quite 
likely that all jets produced by black holes in binary star systems |0] or in AGN's |Q 
Q are initially relativistic. 

Recent jet acceleration and collimation schemes require dynamically strong 
magnetic flelds close to the central engine. Numerical studies ^^ |1^ ||l^ [Q 
show that the jets created in this fashion pass through slow magnetosonic, Alfvenic, 
and fast magnetosonic critical points. The ultimate jet velocity may depend on the 
conflguration of the magnetic fleld Q, and the jets accelerate up to asymptotic speeds 
that may be only a few times the Alfven speed at the Alfven point - at the Alfven point 
the jet speed equals the Alfven speed [n5|. This basic acceleration and collimation 
process may be the same for all classes of objects that emit jets |jl9 . 

Jets are susceptible to the Kelvin-Helmholtz instability. The K-H instability of 
three-dimensional (3D) jets with purely poloidal or purely toroidal magnetic flelds 



|3 11 0' ^^^ of jets containing force-free helical magnetic fields ^ @ 
has been extensively investigated in the supermagnetosonic limit. This instability 
can lead to the development of highly organized helical structures and ultimately to 
disruption of highly coUimated flow as mass is entrained into the jet from the external 
environment 123]. In general, spatial or temporal growth rates associated with the K-H 
instability are found to increase as the magnetosonic Mach number decreases provided 
the jet is super- Alfvenic. Unlike purely fluid flows which are unstable when subsonic, 
the poloidally magnetized jet is predicted to be nearly completely stabilized to the 
K-H instability when the jet is sub- Alfvenic [p3. Here in § 2 and § 3 we consider the 
types of organized structure and stability properties that exist on supermagnetosonic 
and trans-Alfvenic non-relativistic jets, respectively. In § 4 and § 5 we consider the 
stabilizing effects of relativistic jet speeds in the trans-Alfvenic and supermagnetosonic 
regimes. 

2. The Supermagnetosonic Non-Relativistic Jet 

Structures that arise on jets can be considered to consist of Fourier components of the 
form 

/i(r, (f), z) = fi{r)exp[i{kz ±n(j)- ut)] 

where flow is along the z-axis, and r is in the radial direction with the flow bounded by 
r = R. k is the longitudinal wavenumber, n is an integer azimuthal wavenumber, for 
n > the wavefronts are at an angle to the flow direction, the angle of the wavevector 
relative to the flow direction is 9 = tan{n/kR), and ±n refers to wave propagation 
in the clockwise and counterclockwise directions around the jet circumference, n = 
0, 1, 2, 3, 4, etc. correspond to pinch, helical, elliptical, triangular, rectangular, etc. 
"normal" mode distortions of the jet, respectively. For a jet with uniform density, 
temperature and axial magnetic fleld (top hat proflle) and a uniform external medium, 
the propagation and growth or damping of the Fourier components is described by 
a dispersion relation p^. When the flow is supermagnetosonic, each normal mode, 
n, contains a growing single "surface" wave and multiple "body" wave solutions that 
satisfy the dispersion relation (cf., Figure 2). The growth rate of normal mode surface 
and body waves has a maximum that is inversely proportional to the magnetosonic 
Mach number at a wavelength that is approximately proportional to the magnetosonic 
Mach number. The body wave modes exist only for wavelengths less than some 
maximum wavelength. 

Figure 1 shows the maximum displacements associated with asymmetric surface 
wave modes at the fastest growing wavelength and of symmetric and asymmetric body 
wave modes at the maximum unstable wavelength p6[ p7| . The body modes show 
a reversal in fluid displacement at "null" surfaces interior to the jet surface. For 
example, the flrst pinch body mode has a null displacement surface (fluid does not 
move radially) at jet center. The second pinch body mode shows an additional null 
displacement surface between jet center and jet surface. The third pinch body mode 
shows two null displacement surfaces between jet center and jet surface. Pressure 
fronts inside the jet become more oblique to the jet axis as the body wave mode 
number increases. The radial displacement associated with an asymmetric "surface" 



IS 



wave mode n > is approximately given by (,r,n{r) ~ irni''"!-^)^ where ^*„ i 
the surface amplitude |^. Inside the jet surface the displacement is rotated in 
azimuthal angle relative to displacement at the jet surface in addition to changes 
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Figure 1. Axial cross sections associated with maximum pinch body mode (Bl, B2, 
33) displacements and transverse cross sections associated with maximum helical 
and elliptical surface (S) and body mode (Bl, B2, B3) displacements. 



in the displacement amplitude. Fluid displacements associated with the asymmetric 
body waves show a much larger azimuthal rotation and have larger amplitudes in the 
jet interior relative to the jet surface than for the surface wave [^. Only the helical 
surface and body waves lead to jet displacement off the initial axis. Jet distortions 
comparable to the amplitudes shown in Figure 1 have been observed in non-relativistic 
3D MHD numerical simulations pq |. 

3. The Trans-Alfvenic Non-Relativistic Jet 

When a jet is sub-Alfvenic the helical and higher order (n > 0) surface modes are 
stable, but rapidly destabilize when the jet becomes super- Alfvenic. The stability and 
behavior of the normal mode solutions can be investigated analytically in the limit 
uj —f 0. In this limit the real part of the pinch mode (n — 0) surface wave solution 
becomes 1241 
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In the equation above Va = {B^/inpj 



U/2 



Uj is the jet sound speed, and a„ 



(a^ + V^Y^'^ is the magnetosonic speed. The imaginary part of the solution is 
vanishingly small in the low frequency limit. These solutions are related to fast (+) 



and slow (— ) magnetosonic waves propagating with (u+) and against (u—) the jet 
flow speed u, but strongly modified by the jet-external medium interface. Numerical 
solution of the dispersion relation reveals that a growing solution is associated with 
the backwards moving (in the jet fluid reference frame) solution related to the slow 
magnetosonic wave and the pinch mode surface wave is unstable on sub-Alfvenic and 
super- Alfvenic jets. 

When a jet is sub-Alfvenic the helical and higher order {n > 0) surface modes are 
stable [pOJ pi] and have an outwards moving purely real solution given by 

to/k K u + Va ■ 

On the super-Alfvenic jet in the limit a; — > all higher order modes {n > 0) have 
surface wave solutions given by 



oj/k « u < 1 ± i 



1/2 ■ 



1+77 I 7?l/2 J 

where the jet to external medium density ratio r/ = Pj / Px- Wave growth corresponds 
to the plus sign. Note that in the dense jet limit, i.e., ry -^ 00, this expression reduces 
to uj/k « u±Va, and thus the surface waves are related to Alfven waves propagating 
with and against the jet flow speed but strongly modified by the jet-external medium 
interface. The unstable growing solution is associated with the backwards moving (in 
the jet fluid reference frame) wave. The surface wave speed in the observer frame, 
{oj/k)iieai, and the growth rate, {uj/k)imag, are strongly dependent on the density 
ratio. Higher jet density is stabilizing. 

Numerical simulations of low density 3D MHD jets |gj| verify that "light" jets 
destabilize rapidly near to the Alfven point. In the simulations, the observed rapid 
development of instability leads to rapid mass entrainment and slowing of the "light" 
jets in the trans- Alfvenic region before the jets can become supermagnetosonic. 

4. Trans- Alfvenic Relativistic Jet Stabilization 

The rapid development of instability at the Alfven point and subsequent mass 
entrainment and slowing of "light" jets can be overcome. To see what is required 
for enhanced stability of the helical and higher order normal modes we consider the 
asymmetric normal mode solutions to the fully relativistic MHD dispersion relation 
in the limit w ^ [|9): 

/i. ^ 7^^ Ji , ■ [1 - (1 + vKVAhuf - ivh'KVA/crV/' 



where rj = Wj/Wx, Va = (B^ /AnWj)^^^, 7 is the Lorentz factor, and the enthalpy 
W = p + (piri-)(-P/c^) where F is the adiabatic index. A positive value for the 
imaginary part of the frequency indicates instability. Note that the Alfven speed 
can be superluminal. When [1 — (l-t-7y)(VA/7u)^ — (?7/7^)(Va/c)'*] < Othejet is stable 
to the asymmetric normal modes, and in the limit r/ — > cxd it can be shown that 

,, u±Vnj , Va 

uj/k w — — -^5- and v^j = 



1 ± (UV^/C^) " (1 + VJ/c2)l/2 



Thus, helical and higher order surface waves are related to Alfven waves and, of course, 
the Alfven waves propagate at less than lightspeed. The unstable solution is associated 
with the backwards moving (in the jet fluid reference frame) wave. When unstable 
the wave speed in the observer frame, {uj / k) n^ai ~ [I'^'nli^+l'^'n)]^^ increases and the 
growth rate is reduced as ^r^''^ increases. Thus, relativistic jets are effectively heavier 
and can pass through the Alfven point with less rapid development of instability. These 
effectively heavier jets should be much less susceptible to slowing by mass entrainment 
P3i. 



5. Supermagnetosonic Relativistic Jet Structure 

5.1. Theory & Axisymmetric Numerical Simulation Compared 

Jets that survive the trans-Alfvenic region without serious mass entrainment and 
slowing will become supermagnetosonic with continued expansion. The amplitude to 
which normal mode structures can develop on supermagnetosonic relativistic jets can 
be deduced from numerical simulations of relativistic jets |Q [g^. In the simulations 
the jets are perturbed by strong pressure waves driven into the jet at approximately 
constant spatial interval by large scale vortices in the surrounding cocoon. The 
cocoon vortices induce oblique pressure disturbances in the jet at approximately the 
relativistic Mach angle. This oblique disturbance is at the appropriate angle to couple 
to the third pinch body mode (B3 in Figure 1). 

Figure 2 shows normal mode solutions to the dispersion relation appropriate to a 
relativistic jet simulation pfl] (Lorentz factor 7 = 5, Mach number jMj « 10.5 where 
Mj = u/aj, and the jet resides in a hot tenuous cocoon where Mc = u/ac ~ 1.5) 
for surface (S) and the first three body waves (Bl, B2, B3) of the pinch, helical and 
elliptical normal modes. These solutions are typical of supermagnetosonic jet flows. 
With the exception of the surface pinch mode all wave modes have a maximum in the 
growth rate (dashed lines). The body wave modes have a longest unstable wavelength 
(determined from the real part of the wavenumber where the imaginary part goes to 
zero) which is slightly shorter than the longest allowed wavelength (determined from 
the wavenumber in the limit cu —^ 0). Figure 2 also contains a comparison between 
pressure, axial, and transverse (radial) velocity fluctuations seen in the numerical 
simulation with theoretically derived oscillations associated with the B3 pinch body 
mode at the longest unstable wavelength. The comparison indicates that a maximum 
theoretical pressure fluctuation, 2 > P/Pq > 0, successfully reproduces the observed 
maximum transverse (radial) velocity oscillations seen in the numerical simulation far 
behind the jet front, i.e., 15 < z/Rq < 20. 

5.2. Predicted Asymmetric Relativistic Jet Structure 

The maximum pressure fluctuation criterion (2 > P/ Pq > 0) found from axisymmetric 
numerical simulations can be used to predict the cross section distortions and 
asymmetric structural features that will appear on 3D supermagnetosonic relativistic 
jets. In Figure 3 (top panels) we plot the maximum cross section distortions for the 
elliptical surface wave mode and first elliptical body wave mode at the maximum 
growth rate and longest unstable wavelength, respectively, where we have used the 
dispersion relation solutions appropriate to the relativistic jet shown in Figure 2, i.e., 
u/c — 0.98 and 7 = 5.0. Comparison between these distortions and those allowed 
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Figure 2. Solutions to the dispersion relation for the highly relativistic numerical 
simulation (top panels). The dotted (dashed) lines give the real (imaginary) part 
of the wavenumber as a function of the angular frequency. Pressure and velocity 
fluctuations seen in the numerical simulation just off the jet axis, and pressure and 
velocity fluctuations associated with the B3 pinch body mode just off the jet axis are 
shown in the bottom panels. 



on non-relativistic jets (see Figure 1) shows a considerable reduction in the allowed 
distortion. In general, we have found that the allowed distortion decreases as the 
relativistic Mach number increases at wavelengths near to the maximum growth rate 
(surface wave) or between the longest unstable wavelength and the maximally growing 
wavelength (body waves). At longer wavelengths the surface wave elliptical distortion 
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Figure 3. Maximum transverse cross section elliptical surface and first body wave 
distortions (top panels) , and the accompanying pressures and "apparent" emissivities 
around the displacement contours (bottom panels). The dashed (heavy solid, dotted) 
lines in the cross sections correspond to the dashed (heavy solid, dotted) line(s) in 
the pressure and apparent emissivity plots. 



to the cross section can be larger than that shown in Figure 3 pjlj |. 

High pressure regions accompanying the distortions along with the Doppler 
boosting factor 

D = {7[1 - (vlc)cosd\Y^ 



and an angle 6 = l/j + 60 where S9 represents the perturbation to the initial 
flow angle, 9 = I/7, with respect to the line of sight, allows us to construct an 
"apparent" synchrotron emissivity P^D^ . For the elliptical distortions shown in Figure 
3, |i56'| < 0.1 radian, and we also have found that \S9\ oc I/7. The curves of pressure 
and "apparent" emissivity shown in Figure 3 are taken around the displacement 
contours where the dashed (dotted) lines are at the jet surface (center), respectively, 
and the heavy lines indicate the innermost contour of the high pressure region (surface 
elliptical mode) or highest pressure region (first elliptical body mode). The pressure 
and apparent emissivity plots indicate that two bright filaments will accompany an 
elliptical distortion. These filaments will wind around the jet as the elliptical distortion 
rotates spatially down the jet and can show brightness asymmetry as a result of 
variation in the Doppler boosting factor induced by a varying flow direction provided 
the observer is located at about the beaming angle 9 « I/7. At larger viewing angles 
this apparent emissivity asymmetry is reduced. 

Filaments produced by the first elliptical body mode will reside within the jet and 
will dominate filaments produced by the surface elliptical wave mode whose apparent 
emissivity falls rapidly inside the jet surface. The helical first body mode (not shown) 
would generate a single bright twisted filament within the jet which would exhibit 
apparent emissivity asymmetry, similar to that illustrated in Figure 3 for the elliptical 
mode filaments, as the filament winds around the jet. The helical surface mode which 
twists the entire jet is found to result in little apparent emissivity variation |3l| . 

6. Conclusion 

Kelvin-Hclmholtz instability induced normal mode pinch structures have been 
observed in axisymmetric non-relativistic and relativistic jet simulations, and normal 
mode helical, elliptical, triangular and rectangular structures have been observed in 
fully 3D non-relativistic jet simulations. In numerical simulations large amplitude 
helical and elliptical normal mode structures have been found to lead to significant 
mass entrainment. This mass entrainment slows "light" jets and prevents them from 
propagating to large distances while remaining highly coUimated. The growth of 
large amplitude normal mode structures can be slowed by high magnetosonic Mach 
numbers. However, magnetic acceleration and collimation schemes produce jets which 
pass through a trans- Alfvenic and trans-magnetosonic region in which the growth rate 
of the normal modes can be very large. Numerical simulations show that "top hat" 
profile "light" non-relativistic jets will not survive this region. The helically twisted 
structures that we observe in numerical simulations are on much shorter spatial length 
scales than those of twisted structures seen at parsec and kiloparsec length scales on 
extragalactic jets. 

If magnetic jet acceleration and collimation schemes are to prove viable for the 
production of observed jets that propagate to distances orders of magnitude larger 
than the location of the Alfven point, they must flow through the transmagnetosonic 
region with sufficient stability. Additional stability both linearly and non-linearly may 
be achieved by different density, temperature, magnetic and velocity profiles, through 
rapid jet expansion [ p2[ or by the embedding of a jet in a surrounding fast wind. In 
particular, a higher jet density relative to the surrounding environment will slow the 
development of instability. The high jet density of non-relativistic protostellar jets 
relative to the ambient and cocoon medium behind the jet bow shock should allow 
the jets to propagate with sufficient stability through a trans-Alfvenic region. For 



galactic relativistic jets or cxtragalactic highly relativistic jets magnetic acceleration 
schemes would require that the Alfven speed at the Alfven point be near to or greater 
than lightspeed in order to produce a relativistic jet. In this case, theory indicates 
that relativistic effects can be stabilizing and lead to a scenario in which relativistic 
galactic and cxtragalactic jets can propagate through the trans- Alfvcnic region with 
sufficient stability. 

Once through the trans-Alfvcnic region, continued jet expansion leads to 
supcrmagnetosonic flow. On the supcrmagnetosonic jet relativistic effects serve to 
reduce the distortion amplitude of normal mode structures, and this reduction in 
distortion amplitude should lead to a reduction in the mass entrainment that slows 
and disrupts jet coUimation. Although distortion amplitudes are reduced the predicted 
pressure fluctuations combined with the Doppler boosting factor suggest that twisted 
filamentary structures should be readily visible on relativistic jets. 
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